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Abstract. A physically consistent approach is considered for defining an external magnetic field 
as needed in computational fluid dynamics problems involving magnetohydrodynamics (MHD). 
The approach results in simple analytical formulae that can be used in numerical studies where 
an inhomogeneous magnetic field influences a liquid metal flow. The resulting magnetic field is 
divergence and curl-free, and contains two components and parameters to vary. As an illustration, 
the following examples are considered: peakwise, stepwise, shelfwise inhomogeneous magnetic 
fields, and the field induced by a solenoid. Finally, the impact of the streamwise magnetic field 
component is shown qualitatively to be significant for rapidly changing fields. 

There are few examples in recent history, when things being obvious to particular specialists, re- 
mained unexploited by people working in conjugated fields. These include, for instance, the Fast Fourier 
Transform (FFT), which was originally used by Gauss in 1805 and then several times rediscovered by 
Lanczos and Daniclson in 1942, and Cooley and Tukey in mid-1960fl Another and more specific ex- 
ample is the so-called Savitzky-Golay smoothing filter. The least square method being the base of the 
filter has been formulated hundred years before experimentalists working with spectra started to use 
it to treat their data. The paper that popularized this method for the experimentalists is one of the 
most widely cited papers in the journal Analytical Chemistr^. 

What is considered in this letter is not as far reaching as the two aforementioned examples, never- 
theless, we believe it will help people studying numerically the flow of liquid metals under the influence 
of an inhomogeneous magnetic fi eld. Also, this le t ter is complementary to the discussion about mag- 
netic field models given earlier in Votvakov et al. 2008j . The cited work dealt with a 3D distribution 



of heterogeneous magnetic fields, which are of importance in the case of a magnetic obstacle, while 
the current letter considers simpler 2D distributions, which are needed, for example, for the proper 
description of the flow through a fringing magnetic field. 

Let us recall shortly the effects of magnetic fields on conducting liquid flows. When the induced 
magnetic field is negligible, the external magnetic field interacts with the moving liqu id and produces the 



Lorentz force which brakes the flow in the direction of motion, see e.g. lDavidsonl [200lj . This phenomenon 



1 Press, W. H.; Flannery, B. P.; Teukolsky, S. A.; and Vetterling, W. T. "Fast Fourier Transform." Ch. 12.2 in Numerical 
Recipes in FORTRAN: The Art of Scientific Computing, 3rd ed. Cambridge, England: Cambridge University Press, Page 
498, 2007; see also http://mathworld.wolfram.com/FastFourierTransform.html 

2 A. Savitzky and Marcel J.E. Golay (1964). Smoothing and Differentiation of Data by Simplified Least Squares Pro- 
cedures. Analytical Chemistry, 36: 1627-1639; see also http://en.wikipedia.org/wiki/Savitzky-Golay_smoothing_filtcr 
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is heavily exploited in many practical applications ( Davidsonl 1999 1), such as ele ctromagnetic stirring, 



electromagnetic brakes, and non-contact flow measurements ( Thess et al. 20061 ]). 



In reality, one always deals with an inhomogeneous magnetic field because creating a strong and 
homogeneous magnetic field for experiment al needs remains a hard pra ctical challenge. Despite this 
fact, starting from the pioneering work of Hartmann fc Lazarus! jl937 | many theoreticians love to 
work mostly with a homogeneous magnetic field. It is clear, that a constant magnetic field is already 
responsible for main phenomena such as the formation of the Hartmann and parallel layers. Nevertheless, 
a constant field could not prod uce the well develo ped M-shaped velocity profile frequently used for the 
electromagnetic brake. Due to Kulikovskii |l968| . who showed that the flow under the influence of a 
strong and slowly varying magnetic field can be subdivided in a core and a boundary layer, people 
started to exploit theoretically inhomogeneous magnetic fields. It was demonstrated that the flow goes 
parallel to characteristic surfaces, but all the calculations, includi ng the numerical ones, were carried ou t 
either by neglecting the s e cond magnetic fiel d component, see, e .g. Sterll 199C ^ Molokov fc Reedl 200361 



Molokov fc Reed|2003atJAlboussierd[2004j . [Ku"mamaru et al] 120041 jKumamaru et aLl|2007|jNi et al 



2007j or by employing an especially curvilinear channel iToddl |l968| to match boundary conditions 



In the case of numerical simulations, there is no particular technical problem preventing the inclusion 
of the second component of the inhomogeneous magnetic field; nevertheless this has not been done. 
A probable explanation for this is the absence of suitable and analytically simple models to define an 
inhomogeneous magnetic field. Therefore, there is the need for convenient formulae that could allow 
us to do so, and the goal of the paper is to provide a simple method that can be used in numerical 
studies so that both components can be varied in a consistent way. At the end of the paper we discuss 
briefly the legitimacy of omitting the streamwise component of the magnetic field, and show that in 
some cases this might have been done due to large aspect ratios. 

It is worth to note that there are examples where correct expressions for a het erogenous magnetic fie ld 



have been used for MHD flow simulations. These are 2D numerical pap ers bv lCuevas et al. [2006a .6] 
who correctly applied formulae taken from the book of McCaie |l977 |. Those 2D simulations hired 



only the transverse component of th e magnetic field while t he sec ond component played no role. More 



2007, 2008], where all three components of the 



sophisticated 3D cases are given in IVotvakov et 
magnetic field are taken into account. 

The necessity to have at least two nonzero components of the inhomogeneous magnetic field (hereafter 
denoted as MF) follows directly from the requirement that, in the flow region, an externally applied field 
must be simultaneously divergence-free V • B — and curl-free V x B = 0. Thus, if the transverse MF 
component varies along the streamwise coordinate, the streamwise MF component must vary consistently 
along the transverse coordinate. (The spanwise component can be neglected without violation of the 
physical correctness, when the magnetic field is two dimensional.) A vector field that is simultaneously 
divergence-free and curl-free is known as a Laplace vector field and it can be defined in terms of the 
gradient of any function r\ which is harmonic in the flow region, B = — V77, Ar\ = 0. This fu nction 77 



is call ed the magnetic scalar potential. Although this is the most general approach, see e.g. McCaid 



19771 ] . it is not quite convenient because boundary conditions for rj must be defined for each specific 
case. However, what we need is a magnetic field which either vanishes or goes monotonically onto a 
constant level far away from the central point where the intensity of the field is maximal. So, among all 
the possible harmonic functions, we may select those that do not vary at far distances, and this forms 
the basis for the proposed methodology. 

A simple, flexible, and physically consistent way to define an inhomogeneous MF for parametric 
numerical needs is based on the magnetic field induced by a single magnetic dipole. This field is local 
in space and its magnetic scalar potential, which below is called " an elementary potential" , belongs to 
harmonic functions. The whole magnetic field is the sum of local fields from the single magnetic dipoles, 
therefore, the whole magnetic field can also be represented as the gradient of a scalar function, i.e. as a 
whole scalar potential. In other words, the single magnetic dipole can be taken as an elementary unit 
in the appropriate spatial distribution of magnetic sources. The field B'(r, r') creat ed at r = (x,y,z) 
by a single dipole m = (0, 0, m) located at r' = (x 1 , y', z') is given bv I Jackson 19991 ] : 



B'(r,r') = V 



V- 



1 
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where we have used few vector identities and omitted the constant /xq/(47t). Let the dipolcs be distributed 
in the finite region ft outside of the flow. Then, it is easy to see that the elementary rf and the whole 
scalar potential r\ are given by: 

V (x,y,z,x ,z ,z) = -m—^y— —J =y— jTp (2) 

f f JYliz — z') 

ri(x,y,z) = / r) (x,y, z,x' , z , z')dx'dy'dz' = / -nr-dx'dy'dz' . (3) 

J a Jn l r — r r 

This is the simplest model of a permanent magnet occupying the region fl with the constant magnetic 
dipolc distribution m. By taking various fl, one can define any desired inhomogeneous MF for parametric 
numerical needs and mimic roughly various magnetic systems. Below in FigQ]wc give few simple examples 
of two-dimensional, i.e. B y (x, y, z) = 0, fields constructed in this way to represent typical magnetic field 
configurations. Everywhere, x (left to right) is strcamwisc, y - spanwise (perpendicular to the plane of 
Fig- HI), and z (down to up) is transverse coordinate. 

Here it is worth noting that any two-dimensional field can be conveniently expressed through a 
complex representation. Introducing the variable C = x + iz, and its conjugate £* = x — iz, one may 
define a complex function 

m = B x +iB z , (4) 

where B x = Re(/3), and B z = Im(/3). The function /3(C) is given below as well. 

Linear chain of equal dipoles oriented parallel to z direction and infinitely extended in the y direction 

at (xq, zq), fi = {x = .To, — oo < y < oo, z = zq\: 



v(x,y,z) 



m(z — z')S(x' — xq)S(z' — zo)dx'dy'dz' 
. -oo ((x-x>y + (y-y>) 2 + (z-z>) 2 ) 3 / 2 
2m(z — zq) 2m cos ip 



(5) 



(x - x ) 2 + (z - Zq) 2 



dr)(x,y,z) Am(x - x )(z - z Q ) 2msin2<^ 

B x (x, y, z) = — = — — —2 = , (6) 

ox ((x - x ) 2 + {z- zq) 2 ) I 

dr](x,y,z) _ 2 m ((z - z ) 2 - (x - x ) 2 ) _ 2mcos2y 
&z\x,y,z)- — 2 - , (() 

oz {{x - x ) z + (z - z ) z ) I 



here, I and tp define the cylindrical coordinate system (x — xq) = lcos(p, (z — zq) = I sirup, and the 
angle ip is taken as counter-clockwise. Figflji shows an example of the magnetic field induced by two 
symmetrically located chains of dipoles in the rectangular channel. 

The sheet of magnetic dipoles, Fig[Tf>, is made by stacking the linear chains described above along 
the z axis, f2 = {x = Xq, — oo < y < oo, zq < z < oo}. Therefore for this case, the magnetic field is 
obtained by an integration of Eqs. (0[8]) over z, from zq to infinity: 

r](x, y, z) = —m In {(x — xq) 2 + (z — zq) 2 ) = 2m In — (9) 
2m(x-x ) 2mcosip 

Bx{x > y > z) = (x - x ) 2 + (z~ z ) 2 = —r-> (10) 

B z (x,y, z) = 2m l;-^ = (11) 

[X - XqY + {Z - ZqY I 

r(/-\ o C - Co 2m 

m = 2m \t^ = W^cfy (12) 

This particular magnetic dipole configuration allows the construction of the final magnetic field through 
a clear geometrical interpretation. Fig.[Tf> shows this construction: the vector Bi is the field created by 
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Figure 1. Inhomogcncous magnetic field created in the duct by linear chains of identical dipoles(a), two identical sheets 
dipoles(6), half space (c), and two magnets (d). Angles ip, ipi, <p2 must be taken as counter-clockwise, m is taken equal 
one. Also, (6) shows a geometrical construction to find the final magnetic field, see Eq. 1101111 



the lower sheet and is directed along the (r — ri), while its length is inverse to h; on the other hand, 
vector B2 is the field created by the upper sheet, it points along the (r — with a length inversely 
proportional to I2', the vector sum B = Bi + B2 is the final magnetic field in the duct. 

The stepwise magnetic field can be obtained by the integration of the sheet fields from the edge 
point x\ up to infinity, (J? = {x\ < x < 00, —00 < y < 00, —00 < z < Zi\J z% < z < 00}), see 
these coordinates in Fig. [TJ:). To exclude the divergence for B x at x — ► 00, we add a second magnetic 
half-space, then: 



B x {x,y,z) = m In- -—— — = 2mln — 

(x - xxY + {z- z 2 r h 



z — Z\ 



~~2 



B z (x, y, z) = 2 m ( arctan + arctan ■ 

\ x — Xl x — X\ 



2mip, 



/3(C) = 2 m In 



C-Ci 

C-C2 



(13) 
(14) 
(15) 



This, so called fringing, magnetic field is a candidate configuration for the proper study of the trans- 
formation of the Poiscuillc velocity profile to the Hartmann one, and the dependence of a transition 
length as a function of Re and Ha. To our knowledge there are no papers to date that apply a physi- 
cally correct fringing magnetic field in numerical simulations, ev e n though this field is one of the most 
heavily studied in liquid metal f l ows, see Molokov fc Reed 2003 b \ , [Alboussier 1 [2003,| Kumamaru et al\ 



2004| . lKumamaru et all [2007| . lNi et all [2007j . Below we derive qualitatively what is the effect of the 



streamwise magnetic component of the fringing field on the transverse pressure distribution in the duct. 
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The shelfwise magnetic field, created by a homogenous dipole distribution, is obtained if fl is restricted 
in the x direction (J? = {x\ < x < x 2 , — oo < y < oo, -co <z<ziUz 2 <z< oo}). Fig. [TJi shows 
this configuration: 



B x (x,y,z) = ^lnjii^i 

m ((a - xQ 2 + (z - zi) 2 )((z - z 2 ) 2 + (z - zQ 2 ) 
2 n ((a: - xx) 2 + (z - z 2 ) 2 )((a; - x 2 f + (z + z 2 ) 2 ) 
B z (x,y,z) = m(ipi + (p 2 ) 

( Z — Z 2 Z — Z 2 Zi — Z Zi — Z 

= m arctan arctan h arctan arctan 

\ X — X\ X — X2 X — X\ X — X2 



(16) 



(17) 



n(/-\ i (C _ Cll)(C — C21) / 10 s 

/3(CHmln (C-C 1 2)(C-C22)' (18) 

By varying the parameters x\ , x 2 , zi , z 2 , one can regulate the width of the central region and the upward 
(outward) gradient of the magnetic field. 

In above examples, the z integration is taken over two half-open regions {— 00 < z < Z\ (J z 2 < 
z < 00}. It is easy to perform the integration up to a finite z- value instead of infinity, however, this 
case is not considered in order to simplify formulae and keep the essence of the method. Moreover, in 
practice, a magnetic system is supported by the yoke made of soft iron. This effectively means a closure 
for the lines of the magnetic field coming from the external side of the magnet (where the internal side 
is adjoined to the duct), therefore, the upper limit of the integration can be safely set as infinity. 

Up to now, we have learned how to define an external magnetic field by means of magnetic dipolcs. 
Usage of magnetic dipole language dictates the following scaling behavior for the magnetic field far from 
the central point: l~ 3 for a single dipole, l~ 2 for a linear chain, l^ 1 for a sheet, and In I for a massive 
body, where I is a distance from the object, (see Fig. QJ. It's possible to make this slope even more 
smooth, i.e. I In I, if, instead of the magnetic domain, we take a cross section of constant current density 
j, taken everywhere equal to one. Physically, it means that instead of using the permanent magnet 
considered before, we shall deal with a conducting cable carrying the electric current of homogenous 
density. The easiest way to understand how we could develop mathematically the cable from the magnet 
is described next. 

First, we may present the magnetic domain as a non-conducting cross-section enabling surface 
currents only. Schematically, it is shown in Fig. [2b by taking six magnetic dipolcs together in close 
contact: all the internal currents from the adjoining magnetic dipoles are mutually compensated because 
they run in opposite directions, therefore, the only non-compensated electric currents are located on 
the bordering surfaces. The next step is to extend the surface current obtained above, in such a way 
that it flows through a finite rectangular region, i.e. through cable of rectangular cross-section. The 
configuration of the rectangular conductor is shown in Fig. [2b., where the corners are (xk, Zi), k, I = 1, 2 
(the diagonal cross in Fig. [2b denotes electric current running in the direction perpendicular the plane 
of the figure). Such a cable can be represented as a part of the top (cables II and IV) and bottom (cables 
I and III) solenoids assembled around the channel, sec Fig. [2b. 

To derive the magnetic field produced by the aforesaid primitive solenoids, we note that the shelfwise 
magnetic system shown in FigQJf, Eqs. (fT6j) - (fT8l) . is equivalent to the solenoids having infinitely thin 
cables located at x = x\ and x — x 2 because all internal currents in this system arc mutually compensated 
in the way explained above by Fig. [5b.. Therefore, if one integrates given by Eq. (fT5|) over the 
complex variable £ = x + iz confined by cable's corners shown in Fig.[2b, then one obtains the following 
complex function for the individual cable A: 



£22 C , c a _ AN , Cu C 

S.21 — S. S,12 S 



P A (0 = (C 2 1 " In ^r-7 + (Cil - In ^— ^, (19) 



where A = I, II, III, IV should be taken for different cables, and = x^ + izf, k, I = 1, 2 are for the 
coordinates of the corners of cable A, see Fig. \Ep. The final complex function for the top and bottom 
solenoids is obtained as a superposition of terms coming from four cables: 

/3(C) = /3 7 (C) + /3 J/ (C) - P IU {0 ~ P IV (0, (20) 
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fU fll flV flV 

fell fe21 fell fejl 




fl fl fill fill 

fell fen fell fell 



Figure 2. Scheme showing an equivalence of the magnetic domain and a body with surface currents (a), notations used 
to write down magnetic field from a cable (b), scheme of two primitive solenoids located under and above duct(c), the 
cables are supposed to encircle at y = ±00, see text and Eq. 1 1191 1201 1 



where terms from cables I and II are positive and for cables III and IV are negative because the direction 
of the electric current in cables I and II is opposite to the direction of current in cables III and IV. 
The final f3((), Eq. [2Ql is the linear superposition of the products of logarithmic and rational functions 
depending on a complex variable £ = x + i z. Then, by means of algebraic transformations, the obtained 
/3(C) should be presented as a sum of its real and imaginary parts (c/. Eq. [3}, so that B x = Rc[/3(£)] 
and B z = Im[/3(£)]. Although they have been obtained, the mathematical expressions for B x and B z 
are quite cumbersome and for the sake of brevity they are not reported here. 

To conclude, we show briefly what is the impact of the streamwise field component on a liquid metal 
flow in the case of the fringing magnetic field. As it is widely accepted, inertia and viscosity vanish in the 
core of a duct MHD flow, therefore, the pressure distribution, p, is governed in the core by the Lorentz 
force, Vp = j x B, where j are the induced electric currents. This means that Vp is perpendicular 
to B, i.e. the pressure contour lines arc matched with magnetic field lines. Thus, if one employs the 
transverse field component only, the pressure contour lines are straightened in the transverse direction, 
while they must be curved in accordance with the curvature of the external magnetic field, see Fig. [3] 
Curvature of the magnetic field lines depends on the aspect ratio d/h, where d (h) is the height of 
the duct (magnetic gap). As d/h decreases, the curvature vanishes, however the inward gradient of 
the magnetic field decreases as well. Hence, in order to study the effects of a rapidly varying fringing 
magnetic field, the streamwise field component must be taken into account necessarily. 

As a paradigmatic exam ple, the wel l -know n nuclear fusion benchmark fringing magnetic field case 



based on the experiment of iReed et all [1987j (Ha = 6569, N = 10824, c = 0.027) has been recently 
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Figure 4. Dimensionless transverse pressure differenc e for the so called Benchmark fringing magnetic field case. Comparison 
is given for exper i ment al results l|Reed et al\ [l987j |. circles), 3D numerics with a non-consistent fringing m agnetic field 
(dashed. lNi et al\ |2007j |1. and 3D numerics with a consistent magnetic field (solid. lAlbets-Chico et al\ [20091 ]) The results 
are clearly improved when the consistent magnetic field is considered. 



studied using a complete numerical resolution of the Navier-Stokes equations bv lNi et al. 2007]. Details 
regarding the shape of the fringing magn etic field, geometry and boundary con ditions for this case are 



cle arly exposed in both the experimental Reed et al. 1987 1 and the numerical iNi et al. 2007 1 works 



Ni et 



20071 ] presented 3D numerical results in very good agreement with the experimental data 
although a non negligible under-prediction (« 16%) of the peak transverse pressure difference was 
also reported. Ni et al. applied a tanh-based fitting function to approximate the main component 
of the experimental mag netic field, while the additional components were neglected. More recently, 
Albets-Chico et al\ [20091 ] have also addressed this case by means of a complete numerical resolution of the 



governing equations (considering a quasi-static approximation for the induced magnetic field) when using 
a nodal-based non-structured finite-volume code. Details regardin g the form of t he addressed Navier- 
Stokes equations and the assumed simplifications can be found in lNi et aD [20071 ] . as both works have 
employed essentially the same methodology. Additionally, Albets-Chico et al. have analyzed the effect of 
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the consistency of the magnetic field while developing a mathematical technique to generate consistent 



magn etic fields from experimental fits. Further, details will be available soon in lAlbets-Chico et al 
2009], 



Fig. d] presents the effect of the consistency of the magnetic field in the transverse pressure difference 
results. When using a consistent magnetic field, the under-prediction for the peak is reduced to approx- 
imately 9%, which clearly demonstrates how the consistency of the magnetic field plays an important 
role related to the transverse pressure diff erence in such a case. It is interesting to note that Albets-Chico 
et al. obtained a perfect agreement with Ni et al. 20071 ] results when using a non-consistent magnetic 



field (based on the same tanh-bascd fitting function as the one used in Ni et al.). Finally. Albets-Chico 
et al. explain the still remaining undcr-prcdiction in terms of the experimental results in the nature 
(slope and order of accuracy) of the fitting function. 
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